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ABSTRACT 

Light  distributions  of  a  plane  wave  refracted  by  a  microlens  are  calculated  using  Kirchhoff  vector  diffraction 
theory.  Numerical  results  for  one  and  two-dimensional  beam  profiles  and  the  onset  and  effects  of  spherical 
aberration  and  coma  are  investigated  for  different  lens  parameters. 

Keywords:  vector  diffraction  theory,  microlenses,  aberration 


1.  INTRODUCTION 

Ray-tracing  and  geometrical  optics  are  traditionally  used  to  locate  the  focal  regions  of  lenses.  Although  these 
techniques  are  useful  for  quick  calculations  of  focal  distances  and  the  magnitudes  and  regions  of  spherical  aber¬ 
ration  and  coma,  they  are  limited  if  the  vector  components  of  the  light  fields  and/or  detailed  beam  distributions 
are  desired.  For  example,  in  the  field  of  optical  trapping  and  manipulation  of  cold  atoms,  it  is  desirable  to  have 
precise  knowledge  of  the  field  components  of  tightly  focused  laser  beams.  The  trapping  of  cold  atoms  using  the 
dipole  radiation  force  was  first  observed  in  1978,  by  Bjorkholm,  et  al..1  Recent  work  by  Dumke  et  al.2  used  a 
microlens  array  to  focus  laser  light  for  the  trapping  of  cold,  neutral  85 Rb  atoms.  In  this  work  we  apply  vector 
diffraction  theory  to  the  special  case  of  a  microlens  in  order  to  calculate  light  distributions  and  beam  profiles 
in  the  focal  region.  Using  the  radiation  and  vector  diffraction  theory  of  Stratton,3  the  diffraction  integrals  are 
numerically  evaluated  using  average  desktop  computers.  The  onset  and  effects  of  spherical  aberration  and  coma 
on  one  and  two-dimensional  beam  profiles  is  investigated  for  various  microlens  and  light  parameters. 


2.  THEORY 

The  vector  electromagnetic  fields  within  a  given  source-free  volume  of  space  can  be  determined  from  the  vector 
analog  of  Green’s  second  identity  using  Stratton  and  Chu’s  method.4  The  fields  at  point  r  can  be  expressed  in 
terms  of  the  field  values  on  a  surface,  S,  surrounding  the  point  r  using  the  formulae  of  Rothwell  and  Cloud5  as 


—icon  yh'  x  H  (f ' 

H(r)  =  - J  da  [zwe(n'  x  E(r '))  G  +  (V  x  H  (f '))  x  V'G  +  (ri  ■  H  {r  '))  V'g]  ,  (2) 

where  n'  is  the  normal  to  the  input  surface  pointing  into  the  volume,  r  '  denotes  a  point  on  the  surface  S  and 
da'  is  the  surface  area  element  pointing  in  the  direction  of  n  ' .  The  integrations  extend  over  the  input  surface  S 
and  the  values  of  E  (r  ')  and  H  (r  ')  in  the  integrands  are  those  evaluated  ‘just  inside’  the  volume  of  space  after 
the  input  surface.  The  Green’s  function,  G,  used  in  Eqs.  1  and  2  is  given  by 


f)G+(h'xE  (r  '))  x  V'G  +(h'-E  (r  '))  V'g]  ,  (1) 


-ikxp 
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(4) 


where  fci  =  27rni/A,  77 1  being  the  refractive  index  of  the  medium  between  r  and  r',  and 


P 2  =  I r(x,  y,z)  —  f '  (x\  y',  z')  | 


When  the  surface  separating  the  two  media  is  spherical  in  shape,  it  is  convenient  to  use  a  spherical  coordinate 
system  to  evaluate  the  field  integrals.  For  an  incident  plane  wave  traveling  parallel  to  the  0-axis,  the  values  of  the 
electromagnetic  held  components  E  (r  ')  and  H  (r  ')  ‘just  outside’  of  the  spherical  surface  can  be  obtained  from 
the  surface  curvature,  refraction  effects,  and  the  appropriate  Fresnel  transmission  coefficients.  If  E0  denotes  the 
amplitude  of  the  electric  held  of  the  incident  plane  wave,  assumed  to  be  polarized  along  the  ^-direction, 

E(r')  =  E0e-ip2COse'  (tp  sin  9't  cos  <j>'f  +  Tv  cos  0't  cos  <P'0  -  Ts  sin  ,  (5) 

and 

H  (f ')  =  E0e~zp2COsB  (Ts  sin  0't  sin  <jtr  +  Ts  cos  0't  sin  <j>'0  +  Tp  cos  </>'<£)  ,  (6) 

where  Z0  =  P2  =  27T772-R/A,  and  R  is  the  radius  of  curvature.  O'  and  cp'  denote  the  angular  coordinates  of 

the  point  r  ' .  Tp  is  the  Fresnel  transmission  function  for  p-polarized  light, 


m  2ni  cos  O' 

j1  = _ - _ 

p  ni  cos  O'  +  7i2  cos  0't  ’ 

and  Ts  is  the  Fresnel  transmission  function  for  s-polarized  light, 

2771  cos  O' 

S  772  COS  O'  +  77 1  COS  6't 

The  angle  of  transmission,  0’t ,  is  given  through  Snell’s  Law,  or 


=  sm 


(7) 

(8) 

(9) 


2.1.  Normally  incident  light 

The  theoretical  setup  for  this  investigation  is  illustrated  in  Fig.  1.  A  plane  wave  is  traveling  within  the  microlens 
medium  having  a  refractive  index  of  ri2  and  incident  upon  the  spherical  exit  surface  of  the  microlens.  The 
curvature  of  the  spherical  exit  surface  is  assumed  to  have  a  maximum  half-angle  of  0max.  Beyond  the  exit 
surface,  the  refractive  index  is  assumed  to  be  77i,  where  ri2  >  rj  i . 

For  incident  plane  waves  traveling  parallel  to  the  z-axis,  the  electric  and  magnetic  fields  can  be  expressed  as 

E(r)  =  AU(r) ,  (10) 


and 


where 


and 


with 


H  (r)  =  A^V  (r) , 

£ O 


U(r) 


sin  O' dQ' d<p>  [«i  +  U2  +  M3] 


V(f) 


sin  O'dQ'dfj)'  [ui  +  V2  +  M3] 


A  = 


-ipiTp 

Air 


(11) 


(12) 


(13) 

(14) 
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where  pi  =  2nniR/\. 

The  angle  Om  is  the  maximum  half-angle  of  the  aperture,  and  ui,  U2,  U3,  Vi,  fa  and  fa  represent  the  three 
vectors  of  the  three  individual  terms  of  E  and  H  in  Eqs.  1  and  2,  respectively.  Each  of  these  vectors  have  the 
following  components: 


^ lx 

= 

faGi  (cos  9'  cos2  </>'  ■ 

■f  ti  cos  0't  sin2  (/>') 

Uly 

= 

faGi  sin  </>' 

cos  4/  (cos  9'  — 

fa  cos  0't) 

(15) 

Ulz 

= 

-fa  Gi  sin' 

1  9'  cos  cf>' . 

Il2x 

= 

/2G2  [(zi  - 

-  4)<Xl y 

-  (2/1 

-y'i)ot\z\ 

U2y 

= 

/2G2  [(xi  ■ 

—  Xi)  Qliz 

-  («i 

-  4)ai4\ 

(16) 

U2z 

= 

faG2  [(2/1  - 

-y'i)ot  ix 

-  (xi 

—  Xi)  Q:iy]  , 

V3x 

=  /2G2 

sin  9[  cos 

<4  (xi 

-x'l) 

u3y 

=  S2G2 

sin  0't  cos 

4  ( 2/i 

-yfa 

(17) 

U3z 

=  /2G2 

sin  9[  cos 

4  fa  1 

-  4) , 

Vlx 

= 

faGi  sin  4/ 

COS  (/>'  (ti 

cos  O' 

—  cos  0't) 

Vly 

= 

faGi  (ti  cos  O'  sin2  4 

/  +  cos  0't  cos2  4) 

(18) 

Viz 

= 

faGi  (ti  sin  O'  sin  ej>') , 

V2x 

= 

/2G2  [{zi  - 

~  4)  Ply 

-  ( yi 

-y'i)Piz] 

V2y 

= 

/2G2  [(xi  ■ 

~  xi  )01Z 

-  fai 

-  4) Pix] 

(19) 

V2z 

= 

/2G2  [(2/1  - 

-  y'faPix 

-  (xi 

—  xl)  Ply\  1 

V3x  =  hG2h  sin  6't  sin  4  (xi  -  xi) 

V3 y  =  /2G2G  sin  0[  sin  4  (j/i  -  y[)  (20) 

V3 z  =  f2G2ti  sin  0[  sin  4  fai  -  z[)  . 


The  function  fa  is  the  phase  of  the  electromagnetic  field  at  the  point  f  ',  or 

fa  =  e-iP*c°s<>\ 

and  Gi  and  G2  are  the  functions: 


Gi  =  4t tRG  = 


and 


Pi 


G2  =  — 


1  dGr  Gi 


=  —  1  +  t 


1 


(21) 


(22) 


(23) 


ip2pi  dpi  pi  \  ip-2pi. 

The  parameter  t\  is  the  ratio  of  the  Fresnel  transmission  function  for  s-polarized  light  to  the  transmission  function 
for  p-polarized  light,  or 

T„  ni  cos  O'  +  n.2  cos  0't 


1  Tp  m  cos  9[  +  n 2  cos  O' 
The  components  of  du  and  /3\  introduced  in  Eqs.  17  and  20  are  given  by: 


(24) 


a\x  =  sin  (// cos  ((/  (ti  cos  6*'  —  cos0't) 

aiy  =  cos  0't  cos2  <j)'  +  cos  9'  sin2  <j>'  (25) 

a\z  =  —t\sin0' sauj)' , 
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and 


Plx 

Ply 
Pi  z 


—  cos  0'  cos2  4> '  —  ti  cos  0[  sin2  </>' 
sin  (f>'  cos  cj)'  ( ti  cos  0[  —  cos  0') 
sin  0' cos  <j/. 


(26) 


The  coordinates  X\ ,  y\.  z\,  x[,  y[,  and  z[  can  be  expressed  in  terms  of  the  spherical  coordinates  as 


Xi 

=  r\  sin  0  cos  4> 

yi 

=  ry  sin  0  sin  cj) 

(27) 

Z 1 

=  n  cos  0, 

x'l 

=  sin  0'  cos  4> 

v'i 

=  sin  0'  sin  < /> 

(28) 

z'i 

=  cosf?', 

where  r i  =  r/R. 


plane  of 
zero  phase 


Figure  1.  Theoretical  setup  for  calculating  light  fields  beyond  the  exit  surface  of  a  microlens.  The  plane  wave  is  incident 
in  the  +z  direction  and  refracted  at  the  spherical  exit  surface  with  maximum  half  angle  9max ,  and  radius  of  curvature  R. 


2.2.  Light  incident  at  an  angle 

Figure  2  illustrates  the  theoretical  setup  for  light  incident  upon  the  spherical  surface  at  an  angle  of  7  with  respect 
to  the  z-axis.  For  ease  of  calculation,  a  second  (double-primed)  coordinate  system  is  introduced  such  that  x" 
and  z"  are  rotated  about  y'  by  the  angle  of  incidence,  7,  with  respect  to  x'  and  z\  and  y"  =  y' . 

All  of  the  equations  described  previously  for  normally  incident  plane  waves  are  valid  for  light  incident  at 
an  angle,  except  for  the  fact  that  they  apply  to  the  double-primed  coordinate  system;  i.e. ,  all  of  the  “primes” 
in  Eqs.  1-29  are  replaced  with  “double-primes.”  The  difficultly  in  using  the  double-primed  coordinate  system 
arises  in  the  limits  of  the  integrations.  The  limits  of  the  integrations  are  determined  by  the  shape  of  the  exit 
surface  in  that  coordinate  system.  In  the  primed  coordinate  system  of  Fig.  2,  the  limits  of  integration  are  0'  =  0 
to  0m,  and  <j>'  =  0  to  27t.  In  the  double-primed  coordinate  system,  the  limits  of  integration  are  a  complicated 
function  of  0"  and  <f>"  due  to  the  fact  that  z"  does  not  pass  through  the  center  of  the  exit  surface  and  the  exit 
surface  is  not  angularly  symmetric  about  r"  for  all  values  of  0" .  Mathematically  and  computationally  it  would 
be  desirable  to  perform  the  integrations  over  the  primed  coordinates.  If  an  expression  is  found  that  transforms 
the  double-primed  coordinates  to  the  primed  coordinates,  then  the  integration  can  be  performed  in  the  primed 
coordinate  system  where  all  of  the  variables  and  functions  of  U  and  V  (Eqs.  12  and  13)  are  expressed  in  the 
primed  coordinate  system. 
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plane  of 
zero  phase 


Figure  2.  Theoretical  setup  for  calculating  light  fields  beyond  the  exit  surface  of  a  microlens.  The  plane  wave  is  incident 
at  an  angle  of  7  with  respect  to  the  +z  direction. 


The  relationships  between  the  primed  and  double-primed  spherical  coordinates  can  be  found  by  first  express¬ 
ing  the  spherical  primed  coordinates  in  terms  of  the  cartesian  primed  coordinates,  or 

x'  =  sin  9'  cos  4>' 
y'  =  sin  9'  sin  <j> ' 

z'  =  cos  9'.  (29) 

Then  the  double-primed  coordinate  system  is  rotated  about  y"  by  an  angle  —7  to  become  the  primed  cartesian 
coordinates: 

x"  =  x'  cos  7  —  z'  sin  7 

n  1 

v  =  v 

z"  =  x'  sin  7  +  z'  cos  7.  (30) 

Substituting  Eq.  29  into  30,  we  get 

x"  =  (sin0'  cos  4>')  cos  7  —  (cos  9')  sin  7 
y"  =  sin  9'  sin  <// 

z"  =  (sin0'  cos  cj>')  sin 7  +  (cos#')  cos 7.  (31) 


Expressing  the  double-primed  spherical  coordinates  in  terms  of  the  double-primed  cartesian  coordinates  results 
in 


and 


9 "  =  cos  1  (cos  9'  cos  7  +  sin  9'  cos  <j>'  sin  7) 


=  tan 


sin  9'  sin  cf>' 


sin  9'  cos  (j)'  cos  7  —  cos  9'  sin  7 


(32) 

(33) 


Equations  1-29  can  now  be  used  for  light  incident  upon  a  microlens  at  an  angle  of  7  if  the  right  hand  side  of 
Eqs.  32  and  33  are  substituted  for  all  of  the  “0'”  and  “</>'”  variables,  respectively,  in  Eqs.  1-29. 


3.  NORMAL  INCIDENCE  CALCULATIONS 

Figure  3  is  a  calculation  of  the  modulus  square  of  the  electric  field  versus  the  on-axis  distance  from  the  center  of 
curvature.  The  parameters  chosen  for  this  calculation  are  the  same  as  those  for  the  microlenses  by  Dumke  et. 
al 2  to  optically  trap  cold  85Rb  atoms.  The  microlens  array  is  an  array  of  spherically  curved  surfaces  on  one  side 
of  a  fused  silica  substrate  (n  =  1.454  at  A  =  780  nm).  The  light  enters  the  other  (flat)  side  of  the  substrate  and 
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exits  each  spherical  surface.  The  spherical  surface  has  a  radius  of  curvature  of  R  =  312.5  ;tm  and  a  maximum 
half-angle  of  9m  =  11.5°.  Using  geometrical  optics,  the  focal  spot  should  be  observed  at  an  axial  distance  of 
3.2  R  from  the  center  of  curvature.  The  peak  intensity  calculated  in  Fig.  3  is  found  at  an  axial  distance  of  3.11 
R.  Figure  4  is  a  plot  of  the  modulus  square  of  the  ^-component  of  the  electric  field  versus  the  radial  position, 
at  the  focal  plane  (z  =  3.11  R).  One  interesting  feature  of  the  on-axis  intensity  is  the  asymmetry  before  and 
after  the  focal  region.  The  on-axis  intensity  modulations  before  the  primary  focal  plane  resemble  the  on-axis 
intensity  characteristics  for  the  diffraction  of  light  by  a  circular  aperture  without  a  focusing  medium.6,7  The 
on-axis  intensity  features  before  the  focal  spot  indicate  that  light  distributions  in  this  region  are  influenced  by 
the  diffractive  nature  of  the  circular  edge  of  the  exit  surface  as  well  as  the  focal  nature  of  the  spherical  shape 
of  the  exit  surface.  To  further  investigate  the  features  before  the  focal  plane,  Figs.  5  and  6  are  plots  of  the 
radial  intensity  at  the  locations  of  the  on-axis  minimum  and  maximum  occurring  just  before  the  focal  plane, 
respectively.  The  beam  profiles  calculated  in  Figs.  5  and  6  very  closely  resemble  those  of  beam  profiles  calculated 
for  diffraction  of  light  by  a  simple  circular  aperture  using  scalar  diffraction  theory,6  and  vector  diffraction 
theory.7,8 


Figure  3.  Modulus  square  of  the  ^-component  of  the  electric  field  versus  axial  position  for  a  microlens  with  2nnR/\  =  3660 
and  a  maximum  half-angle  of  11.5°. 


Figure  4.  Modulus  square  of  the  s-component  of  the  electric  field  versus  radial  position  for  a  microlens  with  2nnR/\  = 
3660  and  a  maximum  half-angle  of  11.5°. 


4.  SPHERICAL  ABERRATION 

Figure  7  is  a  plot  of  the  on-axis  intensity  for  a  lens  with  a  focal  parameter  of  p2  =  27rn2i?/A  =  5  x  105  for  a  variety 
maximum  half-angles.  A  lens  parameter  of  p2  =  5  x  10s  would  correspond  to  a  radius  of  curvature  of  2.8  cm  for 
532-nm  light.  For  p2  =  5  x  105  the  onset  of  the  effects  of  spherical  aberration  are  evident  for  even  a  maximum 
half-angle  of  9m  =  7.5°.  As  9m  is  increased  from  5  to  7.5°  the  maximum  value  of  the  on-axis  intensity  increases. 
Beyond  a  half-angle  of  7.5°,  the  maximum  on-axis  intensity  no  longer  increases  in  magnitude  with  increasing  6m. 
For  9m  =  10°,  significant  deviations  exist  in  the  on-axis  intensity  profile  from  that  of  an  aberration- free  lens  such 
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Figure  5.  Modulus  square  of  the  z-component  of  the  electric  field  versus  radial  position  for  a  microlens  with  2nnR/\  = 
3660,  a  maximum  half-angle  of  11.5°,  and  an  axial  position  of  the  on-axis  minimum  closest  to  the  focal  plane,  or  z  =  2.69  R. 


Figure  6.  Modulus  square  of  the  a-component  of  the  electric  field  versus  radial  position  for  a  microlens  with  2nnR/\  = 
3660,  a  maximum  half-angle  of  11.5°,  and  an  axial  position  of  the  first  localized  on-axis  maximum  closest  to  the  focal 
plane,  orz  =  2.54 R. 


as  that  of  the  9m  =  5°  on-axis  intensity  plot.  Note  that  this  was  not  the  case  for  a  microlens  with  p2  =  3660  and 
Qm  =  11.5°,  discussed  in  the  previous  section,  illustrating  that  the  onset  of  spherical  aberration  is  dependent 
upon  the  ratio  of  the  curvature  of  radius  to  the  wavelength  of  light.  The  higher  the  ratio,  the  earlier  the  onset 
of  spherical  aberration  for  increasing  maximum  half-angles.  In  Fig.  7(b)  the  maximum  half-angle  is  increased 
even  further  to  investigate  the  effects  of  severe  spherical  aberration.  Once  the  maximum  on-axis  intensity  has 
saturated  (around  7.5°  for  p2  =  5  x  105),  the  region  of  higher  intensities  elongates  along  the  axis  of  propagation 
towards  the  exit  surface  of  the  lens.  The  growth  of  the  “focal  spot”  towards  the  exit  surface  of  the  lens  can  be 
explained  using  geometrical  optics.  According  to  geometrical  optics,  if  a  series  of  parallel  rays  are  refracted  by  a 
spherical  surface  then  the  paraxial  rays  focus  at  the  expected  focal  spot.  Noil-paraxial  rays  will  be  refracted  by 
the  spherical  surface  such  that  they  will  cross  the  optical  axis  at  a  point  closer  to  the  exit  surface.  The  farther 
the  initial  ray  is  located  from  the  optical  axis,  the  closer  to  the  surface  of  the  lens  the  ray  will  cross  the  optical 
axis. 

Figures  8(a)  and  (b)  are  plots  of  the  radial  intensity  profile  at  the  focal  plane  of  a  plane  wave  refracted  by  a 
spherical  exit  surface  as  a  function  of  the  maximum  half-angle  on  a  linear  and  a  logarithmic  scale,  respectively. 
The  intensity  axis  of  Fig.  8  is  normalized  to  the  incident  intensity  before  the  spherical  surface.  Once  the  on-axis 
maximum  intensity  saturation  point  is  reached  for  increasing  0m,  no  significant  changes  in  the  radial  beam  profile 
is  observed  on  a  linear  scale  for  increasing  9m.  Figure  8(b)  is  the  same  plot  as  Fig.  8(a),  except  that  the  intensity 
scale  is  logarithmic.  On  a  logarithmic  scale  it  is  apparent  that  the  light  collected  by  the  outer  rings  of  the 
spherical  surface  crosses  the  “focal”  plane  at  off-axis  locations;  i.e.,  the  energy  collected  by  the  outer  regions  of 
the  spherical  surface  is  propagated  to  off-axis  locations  in  the  focal  plane.  Away  from  the  central  focal  region, 
the  intensity  of  a  spherically  aberrated  beam  can  be  orders  of  magnitude  higher  than  that  of  a  non-aberrated 
beam. 
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Figure  7.  Calculated  values  of  the  intensity  (normalized  to  the  incident  intensity)  along  the  axis  of  the  lens  (normalized 
to  the  radius  of  curvature  of  the  lens)  for  various  values  of  the  maximum  half-angle  (9m  (5°,  7.5°,  10°,  15°,  20°,  and  25°) 
of  the  lens  surface  for  a  normally  incident  plane  wave  with  2nnR/\  =  5  x  105. 


Figure  8.  Calculated  values  of  the  focal  plane  intensity  distributions  (normalized  to  the  incident  intensity)  along  the 
axaxis  (normalized  to  the  radius  of  curvature  of  the  lens)  on  a  linear  and  long  scale,  respectively,  for  various  values  of  the 
maximum  half-angle  9m  (5°,  8°,  10°,  15°,  and  25°)  of  the  lens  surface.  The  light  is  a  normally  incident  plane  wave  with 
2-KnR/X  =  5  x  10s,  and  the  ‘focal  plane’  occurs  for  z/R  =  2.984. 
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5.  OFF-AXIS  FOCUSING  AND  COMA 


Figure  9  is  a  calculation  of  the  intensity  of  the  light  field  versus  x  and  y  in  the  focal  plane  (z/R  =  3.15  for 
27T712-R/A  =  104)  and  incident  angles,  7,  of  0,  7  and  15°.  With  an  incident  angle  of  7°  the  effects  of  coma  are 
observed  as  subtle  partial  rings  on  the  far  side  of  the  focal  spot,  and  a  slight  elongation  of  the  focal  spot  in  the 
direction  of  the  incident  angle.  For  an  incident  angle  of  15°,  severe  coma  is  observed  as  a  broadening  of  the  focal 
spot  in  both  the  x  and  y-directions  in  addition  to  conical  beam  fanning  as  expected  according  to  geometrical 
optics.  However,  contrary  to  geometrical  optics  (which  would  predict  a  smooth  intensity  decay  throughout  the 
beam  fanning)  intensity  calculations  presented  here  predict  the  interference  patterns  arising  from  the  wave  nature 
of  light.  The  degree  of  coma  observed  is  also  a  function  of  the  radius  of  curvature  to  wavelength  ratio.  Similar  to 
Fig.  9,  Fig.  10  is  a  series  of  intensity  versus  x  and  y  images,  where  the  only  difference  is  that  the  value  of  p2  has 
been  decreased  from  104  to  3660.  For  smaller  p2,  the  effects  of  coma  for  7  =  7°  have  essentially  been  eliminated, 
and  the  coma  effects  for  7  =  15°  have  been  significantly  reduced,  as  compared  to  Fig.  9. 

Figure  11  is  a  calculation  of  the  beam  intensity  profile  (along  the  x-axis)  for  the  same  parameters  as  that  of 
Fig.  9,  but  a  series  of  angles  of  incident  leading  up  to  7°.  The  relative  ranges  for  each  graph  are  kept  constant  to 
aid  in  the  visual  comparison  of  the  absolute  intensities  and  radial  beam  fanning.  Even  for  an  angle  of  incidence 
of  3°  asymmetries  in  the  beam  profile  are  observed  for  2nnR/\  =  104.  As  the  angle  of  incidence  is  increased,  the 
overall  maximum  intensity  decreases  and  the  width  of  the  focal  spot  increases  as  well  as  the  amount  of  energy 
observed  outside  of  the  primary  focal  spot. 

15 

S  10 

o 

.9  5 

E 

c  0 

o 

1  -5 

§-  -10 

X 

-15 


Figure  9.  Calculated  intensity  distributions  in  the  focal  plane  for  ZirnR/X  =  104  with  9m  =  11.54°,  n  =  1.454,  and 
incident  angles  of  7  =  0,  7°  and  15°,  respectively. 
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Figure  10.  Calculated  intensity  distributions  in  the  focal  plane  for  2nnR/X  =  3660  with  9m  =  11.54°,  n  =  1.454,  and 
incident  angles  of  7  =  0,  7°  and  15°,  respectively. 
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Figure  11.  Effects  of  coma  on  the  radial  beam  profile  of  a  plane  wave  incident  on  a  microlens  spherical  exit  surface  with 
2nnR/\  =  104  and  z/R  =  3.15  at  an  angle  of  (a)  3°,  (b)  5°,  and  (c)  7°.  All  horizontal  and  vertical  scales  are  over  the 
same  ranges. 
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6.  CONCLUSIONS 


A  computational  method  has  been  employed  to  determine  all  of  the  vector  electromagnetic  field  components  for 
a  region  of  space  beyond  a  spherical  microlens  using  vector  diffraction  theory.  Using  this  method  the  held  and 
intensity  distributions  are  obtained  and  the  onset  and  effects  of  spherical  aberration  and  coma  can  be  studied. 
The  onset  and  severity  of  spherical  aberration  depend  on  both  the  maximum  collection  angle  of  the  lens  and  the 
ratio  of  the  radius  of  curvature  to  the  wavelength  of  light.  The  onset  and  severity  of  coma  depend  on  both  the 
angle  of  incidence  and  the  ratio  of  the  radius  of  curvature  of  the  lens  to  the  wavelength  of  the  incident  light. 
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